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ABSTRACT 


The  temperature  fields  are  deduced  for  an  infinite  metal  slab  of 
finite  thickness  produced  by  a  heat  flux  input  which  depends  on  the  radial 
coordinate.  The  formal  solution  is  obtained  by  means  of  transform 
techniques,  and  several  alternative  approaches  to  their  evaluation  are 
presented.  Two  heat  flux  inputs  are  considered  in  detail.  They  are 
the  Gaussian  input  defined  in  Eq.  (E-l),  and  the  uniform  input  over  a 
circular  area  of  the  front  surface.  Such  inputs  might  be  used,  for 
example,  as  models  for  a  study  of  the  thermal-elastic  effects  of  a  laser 
beam  focused  on  a  metal  slab. 

Several  asymptotic  expansion  are  presented  for  the  small-time 
response  of  the  temperature  fields  under  various  assumptions  on  the 
heat  flux  input. 

An  extension  of  the  results  of  Oosterkamp3  has  been  derived  as 
well  as  a  verification  of  the  validity  of  the  approximations  through  the 
derivation  of  an  upper  bound  on  the  error. 
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THE  THERMAL  RESPONSE  OF  A  METAL 
SLAB  TO  A  CLASS  OF  RADIALLY 
DEPENDENT  HEAT  INPUTS 


A .  Introduction 

The  temperature  fields  produced  by  a  heat  flux  input  dependent 
on  the  radial  coordinate  only  are  deduced  for  an  infinite  metal  slab  of 
finite  thickness. 

The  formal  solution  is  obtained  by  the  use  of  transform  techniques, 
and  several  alternative  approaches  to  their  evaluation  are  presented. 

Attention  is  focused  on  a  class  of  input  flux  densities  which  are 
concentrated  at  the  origin,  such  as  would  occur  if  a  laser  beam  were 
focused  to  a  small  spot  on  one  surface  of  the  slab.  Two  examples  of 
this  type  which  have  been  considered  in  detail  are  the  Gaussian  and 
uniform  inputs  defined  in  Eqs.  (E-l)  and  (F-l),  respectively.  Unfortu¬ 
nately,  extensive  numerical  information  can  be  obtained  only  at  the  cen¬ 
ter  of  the  front  surface.  In  terms  of  the  geometry  described  in  Fig.  1, 
this  is  at  p  =  £  =  0.  This  is,  however,  a  very  interesting  point  since 
the  temperature  maximum  for  the  class  of  concentrated  inputs  occurs 
here.  The  temperatures  at  this  point  also  illustrate  some  of  the  dis¬ 
tinctive  features  of  the  problem. 

One  of  the  basic  results  shows  that  for  time  sufficiently  small, 
the  response  is  essentially  a  local  phenomenon.  Specifically,  this  means 
that  if  q(p)  is  a  sufficiently  smooth  function  of  p,  the  primary  diffusion 
is  into  the  slab  while  the  effect  of  radial  diffusion  is  secondary  for  time 
small.  This  in  turn  means  that  the  temperature  response  under  these 
conditions  should  be  given  to  a  good  approximation  by  the  product  of  the 
local  input  multiplied  by  the  temperature  response  for  the  flat  plate  with 
a  uniform  input  on  the  front  surface.  In  this  paper  several  results  of 
this  type  will  be  given,  each  with  a  different  type  of  "smoothness " 
requirement.  Remainders  will  indicate  their  validity. 

The  paper  can  be  divided  into  two  parts.  The  first  part  deals 
with  the  general  solution  to  the  problem,  and  its  various  formulations. 

The  second  part  deals  with  two  specific  inputs  for  which  more  explicit 
results  can  be  obtained.  These  results  are  presented  in  analytical  and 
graphical  form  when  possible. 
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Fig.  1.  This  figure  illustrates  the  geometry  of  the  problem. 
The  heat  flux  input  occurs  on  the  front  surface 
(Z  =  0)  and  depends  on  the  radial  coordinate  only. 
The  rear  surface  (Z  =  i)  is  insulated. 


This  problem  has,  of  course,  been  studied  by  others,  but  the 
main  results  of  the  other  authors  have  been  for  the  semi-infinite  body 
with  a  uniform  input  over  a  portion  of  the  surface  and  for  final  temper¬ 
atures  in  the  cases  where  a  linear  radiation  law  is  included. 

CarsJaw  and  Jaeger1  give  several  references  to  this  type  of 
problem,  and  other  than  these  references  it  has  not  been  possible  to 
find  any  others.  Lowan2  considers  the  semi-infinite  body  with  a  more 
general  time  dependence  than  that  considered  here.  In  the  present 
report  the  input  flux  has  a  step  dependence  with  respect  to  time  and  is 
a  basis,  through  Duhamel's  theorem,  for  the  study  of  the  more  general 
time  dependence  problem. 
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Oosfcerkamp3  considers  a  uniform  spot  input  on  the  front  surface  of 
a  semi  “infinite  body.  He  reasoned  correctly  that  the  temperature  re¬ 
sponse  for  the  semi-infinite  body  would  furnish  first-order  approximation 
for  the  finite  slab  for  time  small.  This  is  because  for  time  small  the 
presence  of  the  rear  wall  is  net  significant,  nor  is  the  effect  of  radial 
diffusion.  He  did  not  present,  however,  an  evaluation  of  this  approxi¬ 
mation,  This  paper  extends  his  result  and  verifies  his  reasoning  by 
giving  a  rigorous  estimate  of  the  error.  While  the  approximation  itself 
does  not  depend  on  the  thickness  of  the  slab,  the  error  estimate 
appropriately  does, 

Thomas4  considers  the  semi-infinite  body  and  the  very  thin  slab. 

His  results  for  the  semi -infinite  body  include  a  linear  radiation  loss. 

The  input  to  his  thin- slab  problem  is  not  a  surface  flux  but  a  constant 
volume  input  restricted  to  a  right~circular  cylinder  domain  normal  to 
the  slab  surface. 

The  motivation  for  the  present  problem  is  similar  to  Oosterkamp's . 
He  was  interested  in  the  thermal  effects  of  electron  bombardment  of  the 
anode  of  an  X-ray  tube,  and  how  to  prolong  the  life  of  the  tube.  In  the 
present  problem  the  interest  is  in  determining  the  effect  of  the  focusing 
of  a  laser5  output  on  a  metal  slab.  The  laser  is  a  device  which  radiates 
an  almost  coherent  electr omagnetic  energy  at  optical  frequencies.  This 
energy  from  the  laser  can  be  focused  to  a  very  small  spot.  A  typical 
diameter  is  one  millimeter.  The  resulting  incident  power  densities  are 
about  1010  watts/meter2  .  The  coupling  of  electromagnetic  energy  to 
metals  is  quite  efficient  at  optical  frequencies,  especially  when  compared 
to  the  coupling  at  microwave  frequencies.  In  some  cases  the  laser  beam 
will  punch  a  small  hole  in  a  thin  metal  slab.  It  is  felt  that  the  first  step 
in  a  study  to  determine  the  cause  of  this  punching  is  an  analysis  of  the 
temperature  fields.  From  this  study,  then,  the  study  of  the  thermally 
induced  stress  fields  can  be  initiated.  6 

30  3 1 

Two  recent  papers  9  *  indicate  that  the  solution  to  this  problem  can 
also  be  of  value  in  improving  the  accuracy  of  the  determination  of  the 
thermal  parameters.  First,  the  effect  of  the  shape  of  the  beam  can  be 
taken  into  account  and  secondly  the  plate  size  can  be  made  such  that 
the  edges  of  the  plate  don't  effect  the  determination. 
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B*  Formal  Solution 


In  this  section  the  formal  solution  to  the  problem  will  be  presented 
in  two  alternate  forms:  a  Green’s  function  form,  and  a  transformed 
Green's  function  form.  Both  are  obtained  by  transformation  techniques 
which  are  quite  natural  for  the  problem.  The  Laplace  transform  is  used 
because  of  the  time  dependence  of  the  input  flux  density;  it  is  zero  for 
time  negative.  The  Hankel  transform  of  zero  order  on  the  radial 
dependence  is  used  because  of  the  cylindrical  symmetry  of  the  input  flux 
density.  These  two  transformations  reduce  the  problem  to  that  of  finding 
the  solution  for  the  plane  slab  problem  with  a  uniform  input  over  the 
entire  surface  and  performing  two  inverse  transforms.  The  use  of  one 
or  the  other  of  the  two  forms  of  the  solution  depends  on  the  input  flux 
density  q(p  )  and  its  Hankel  transform,  Q(\  ).  The  geometry  of  this 
problem  is  shown  in  Fig,  1  and  the  differential  equations  and  boundary 
conditions  are  given  in  Eq.  (B~l): 


(B-l) 


I  _JL  f  9AU  \  a2A'U 

p  9p  i^p  9P  )  +  9£2 


9AU 

8t 


:0, 


0  <  p  <  oo 

°  <  i<  i 

0  <  T 


9AU(p  ,  0,  t  )  _  _  q(p); 

9£ 


0  <  p<  00, 


T  >  0 


=  0 


T  <  0 


9AU  (p  ,  1,  t  )  _  q 

9£ 


0  <  p  ,  T  <00 


0  <  p  <  oo 

AU(p,£,0)  =  0 

0  <  .£<  1 

where  AU(p  ,  £,  t  )  ~  T{r,  z,  t)  •  T0  and  ^  ~  z/ 1 ,  p  =  r/l ;  I  =  slab  thickness; 
t  -  Kt/£2  ;  K  is  thermal  diffusivity,  t  is  time;  q(p  )  -q(r)l/k;  k  is  thermal 
conductivity. 

The  normalized  heat  flux  input  density  q(p  )  will  be  called  the  heat 
flux  input  density  or  heat  flux  input.  The  units  of  the  normalized  heat 
flux  input  density  q(p  )  in  the  MKS  system  are  in  degrees  Kelvin.  Later 
in  the  paper  the  symbol  Q  with  no  argument  will  be  used  to  denote  the 
total  heat  input.  In  this  case  it  will  be  given  by 
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00 


Q  =  2tt 


J  q(r)  r 

p=0 


and  measured  in  watts;  that  is,  Q  is  the  actual  total  heat  input. 


The  formal  solution  using  the  Laplace  transform  and  Hankel 
transform  techniques  is  given  by  the  expression 


o+  +  ioo 


(B-2) 


AUIp'^t)=^T  I 


QOO 


COGil 


S=0^  -  loo  X  -0 


{vjx2  +  s  (1-9} 


.  ebT  .  J0(Kp  )K 


J  X.2  +  s  s inh  \J  \2  <¥  s 

where  Q(X  )  is  the  Hankel  transform  of  the  input  flux  density,  q(p  ),  and 

00 

(B-3)  Q(\)  =  J  q(p')  J0(Xp’)  p1  . 

p  r~o 


From  this  point  there  are  indeed  many  paths  to  follow.  As  can  be  seen, 
there  are  three  integrations  to  perform  and  thus  there  are  a  multitude 
of  possible  approximations  which  can  be  made.  It  is  impossible  to 
discuss  them  all.  Therefore,  the  remainder  of  this  report  discusses  the 
methods  most  advantageous  to  the  needs  of  the  laser  problem  mentioned 
earlier . 


First,  the  integration  with  respect  to  s  will  be  carried  out.  There 
seems  to  be  no  alternative  other  than  this.  Notice  that  if  X  were  zero  the 
bracketed  term  given  below  (from  the  integrand  in  Eq.  (B-2)), 


(B-4) 


H(£,X  ,  s) 


1  cosh  { >jx 2  +  s  (1  -  4)} 
s  \j\2  +  s  sir.h  \|  X 2  +  s 
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would  be  the  Laplace  transform  of  the  slab  problem  for  a  uniform  heat 
flux  input  over  the  entire  front  surface.7  In  terms  of  the  present 
problem  the  input  flux  q(p  )  would  be  a  constant  for  all  p  and  although 
Q(X  )  would  not  exist  in  the  ordinary  sense,  it  would  be  a  functional, 
namely,  6  {X  ),  When  this  is  substituted  in  Eq.  (B~2)  the  result  is  the 
constant  flux  solution.  This  is  simply  a  statement  of  Hankel's  integral 
formula.*  For  X  f-  0  there  is  an  inter esting  relationship  between  this 
problem  and  the  thin-* wire  problem  with  linear  radiation  from  the  lateral 
surface  as  discussed  by  Sommerfeld.  9  Basically,  the  two  problems  are 
quite  similar.  'This  fact  is  apparent  when  one  notes  that  the  differential 
equation  satisfied- by  the  inverse  Laplace  transform  of  H(£,  X  ,  s)  is  the 
same  as  the  thin-wire  problem  heated  uniformly  at  one  end.  The 
variable  X  corresponds  to  the  radiation  loss  parameter  of  Sommerfeld's 
problem.  The  solution  to  the  present  problem  is  then  seen  as  an 
integration  over  all  possible  radiation  loss  parameters  multiplied  by  an 
appropriate  weight  factor.  The  main  difference  is  that  in  the  present 
problem  the  radial  diffusion  is  caused  by  the  unequal  heat  flux  input  on 
the  front  surface  rather  than  a  radiation  loss  from  the  lateral  surface. 
Also  when  no  radiation  loss  is  present  (X  ~  0)  one  expects  to  obtain  the 
flat-plate  solution  with  no  radiation  loss.  This  is  what  is  obtained,  as 
mentioned  earlier. 


Performing  the  inverse  Laplace  transform  one  arrives  at  the 
following  expression  for  the  temperature  field;  this  will  be  called  the 
transformed  Green's  function  representations. 


(B--5) 


where 


AU  ( p  ,  T  )  -  y  Q(X  )  F(  p ,  t  ;  X  )  X 
X  =0 


F(o,e.  t  ;  X)  ^  JG(Xp)  f 

r  ’=0 


and  03 (£/2,  iirT  ')  is  a  theta  function.  10  The  two  representations  of 
the  theta  function  are  given  in  Eqs.  (Boa)  and  (B-6b). 

°°  2  2 

(B-6a)  0 3  ^  ,  in  ij  -  1  +  2  ^  cos(mr£)e  n  77  T 

n  =  l 


1  107-14 


6 


(B- 6b) 


Again  it  is  seen  that  for  X  =  0  the  t  1  integration  yields  the  temperature 
fields  for  the  constant  heat  flux  input  over  the  entire  front  surface. 

For  X  $  0  the  integration  with  respect  to  t  1  can  be  performed,  but 
it  is  not  particularly  advantageous  to  do  so  at  present.  If  the  integration 
is  carried  out,  one  of  the  resulting  series  is  independent  of  t  .  It  can  be 
summed,  11  but  because  of  its  complexity  it  is  not  particularly  useful. 


Notice  that  in  Eq.  (B~5)  there  are  still  three  integrations  to 
perform.  The  improvement  of  Eq,  (B-5)  over  Eq,  (B-2)  is  that  the 
theta  function  has  two  representations.  One  of  these  converges  well 
for  t  small,  Eq.  (B-6b),  and  the  other  for  r  large,  Eq.  (B-6a). 


Equation  (B~5)  can  be  rewritten  in  a  form  calling  direct  attention 
to  the  fact  that  there  are  three  integrations  to  perform.  This  form 
displays  the  Green’s  function  for  the  problem,  and  for  some  purposes 
it  is  more  useful  than  that  given  in  Eq,  (B-5).  This  form  is  given  by 


(B-7)  AU(p  ,  £>  t  )  =  f  q(p  )  G(p  ,  p  '  ,  %,  t  )  p 

P  ~o 


where  the  Green's  function  is  given  by 


(B-8) 


T  00 

G(p  ,  p  §,  t  )  =  J  ^  J0(Xp  )  Jo<xP  ')  ®3 

t'=0  \  =0 


X. 


Thus  Eqs,(B-5)  and  (B-7)  give  two  alternate  representations  of  the 
solution.  The  use  of  one  or  the  other  depends  on  whether  the  input 
flux  density  q(p  )  or  its  transform  Q(X  )  is  more  convenient  to  use. 

There  are  several  ways  to  deal  with  G(p  ,  p  £,  t  ),  and  apparently 
only  one  way  to  deal  with  F(p  ,  £,  t  ;  X  )  without  explicit  information  about 
Q(X  )  or  q(p).  Unfortunately  not  all  of  the  results  obtained  for  the  Green’s 
function  G(p,  p 1 ,  r  )  and  its  transform  F(p,  t  ;X  )  are  useful  on  the 
special  cases.  This  is  because  most  of  the  significant  results  for  the 
Gaussian  and  uniform  inputs  depend  explicitly  and  strongly  on  the 
particular  form  of  the  input  flux  density  or  its  transform.  For  each  form 
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F(p,  t  ;X  )  and  G(p,  p  1  ,  §,  t  )  it  is  desirable  to  have  both  an  expansion 
which  converges  well  for  time  small  and  one  which  converges  well  for 
time  large.  These  results  are  given  for  F(p,  £,  t  ;X  )  in  Section  C  and 
for  G(  p,  p 1  ,  r  )  in  Section  D. 


C.  The  Transformed  Green's  Function 


Two  general  results  based  on  Eq.  (B-5)  are  deduced  in  this  section. 
One  is  an  expansion  of  F(p,  £,  t  ;X  )  for  t  large  which  is  the  basis  for 
appr oximations  to  AU(p,£,t)  for  time  large.  The  second  is  an  asymptotic 
expansion  for  time  small,  with  some  restrictions  on  the  input  flux  density, 

q(p). 


The  large-time  approximation  will  be  deduced  first.  The  first  step 
is  to  perform  the  t  1  integration.  This  yields  the  following: 


(C-1) 


AU(p,£.T)  =  J  Q(X) 
X  =0 


1-e 


2  _  ! 


“X  “  T 


Y"  cos  mr  £ 

+  2  )  — - — -  J0 

n4  ^  +“***] 


3(Xp)X 


-I  Q(x>  2 1 

X  =0  _  n=l 


cos  mr£ 
X2  +  n2  tt 2 


exp  -(X2  +  n2  ir  2)i 


JQ(Xp  )X 


The  first  observation  is  that  there  are  three  integrations  left.  These 
are  as  follows: 


(i) 


00 

f  Q(X) 

1-e"^ 

J 

^  =0 

X2 

J0(X  P)X 


(ii) 


X  =0 


X  +  nc  ir 


2  „  2 


(iii) 


s 


X  =0 


Q(X)A2±HzL  jo(V|x 


X2  +  n2  TT2 
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To  ask  that  all  of  these  be  expressable  in  a  simple  fashion  for  some 
Q(X )  is  too  much.  The  second  observation  is  that  as  r  becomes  larger, 
the  last  integration  becomes  small  quite  rapidly.  It  is  easy  to  verify 
that  this  last  integration  in  Eq.  (C-l)  is  bounded  by  the  quantity 

-7T2  T 

(C-2)  R  <  -  ,  0<p<«5,  0<  e<  1  . 

This  expression  does  indeed  vanish  rapidly  with  increasing  t,  and 
Q(o)  is  proportional  to  the  total  heat  flux  input  in  watts  to  the  slab, 


Q(o) 


Q 

2ir£k  , 


where  Q  is  the  total  heat  flux  input.  In  conclusion,  then,  the  temperature 
change  for  large  t  is  given  approximately  by  the  expression 


(C  -  3) 


AU(p, 


r°° 

r  \  2  _  oo 

1  “  \  T  r -  £ 

1-e  \  cos  nTT  9 

) Q(x> 

x2  +2  2  X2  +  n2  it2 

=  0 

L  n=l  J 

J0(Xp)X  +Rj( 


where 


R  < 
1 


Q_ 

fk 


e 

12tt  t 


and  where  Q  is  the  total  heat  flux  input.  This  is  the  "large -time"  approxi¬ 
mation  based  on  the  F(|S),£,t;\)  representation.  It  yields  a  term  which 
depends  on  t  and  a  term  which  depends  only  on  the  space  coordinates. 

It  is  assumed  that  the  remaining  integration  can  be  performed  in  a  form 
suitable  for  numerical  computation. 

The  first  integration  (i)  presents  some  difficulties  for  all  values 
of  t.  The  problem  arises  in  the  treatmeatment  of  the  indeterminate  form 


X2 


at  X  =  0.  It  happens  that  in  this  problem  Q(X  )  is  an  even  function  of  X 
(see  Eq.  (B-3)),  and  further  Q(X  )  at  X  =  0  is  not  zero.  Thus  if  one 
tries  to  break  this  integral  up  into  the  sum  of  two, 
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( C  -4) 


oo 


r- 

'  “X2  T  " 

1  -  e 

X2 

c  =0 

- 

JQ(Xp)X 


X  =  0 


JQ(Xp  )X 


-r 


Q(M  — t 


it  is  seen  that  individually  the  integrals  do  not  exist.  Thus  the  only  way 
it  seems  to  be  possible  to  handle  this  integration  is  by  the  device  used  to 
obtain  this  representation, 


(C-5) 


00 

Q(X) 

J 

X2 

x=o 

J0(Xp)X 


■  r  r  <»> 

t'=0  X  =  0 


e 


-XV 


JQ(Xp  )X  . 


This  problem  will  be  met  again  when  the  Gaussian  input  is  studied  and 
will  be  handled  in  this  manner.  Later  the  expression  in  Eq.  (C-3)  will 
be  used  to  deduce  some  information  on  the  temperature  fields  for  the 
uniform  and  Gaussian  inputs. 

Next  an  asymptotic  expansion  for  the  temperature  field,  £U(p,  £,t), 
as  t  -►  0  will  be  presented.  This  expansion  is  discussed  more  fully  in  an 
earlier  report.12  The  sufficient  assumptions  for  the  expansion  are  as 
follows : 

d2N 

(C-6)  (a)  - — •  q(p  ) 

dp*2 


exist  for  some  N,  and 


<»  I" 

X  =  0 


|Q(x)|x 


2N+2 


exists . 


These  two  conditions  constitute  the  smoothness  requirements  for  this 
development. 


Under  these  assumptions  the  X  integration  in  Eq.  (B-5)  can  be 
performed  to  yield 


(C-7) 


where 
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/ 

ITT  T 


I(p.  t') 
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l{p>T  ) 


q  ( p )  +  ^ 


n  f  1  d  d  \ n 


n !  p  dp  dp 


q(p)  +  Ra 


■p  ^  - 

(N+l) ! 


|Q(M|  X 


It  is  noted  that  Ra  is  o(t  1NJ)  as  t  “*0.  The  final  result;  an  asymptotic 
expansion  to  N  terms  as  t  ^0  of  ^U(p;£>t),  is  given  by  the  expression 

N  x 

(C-8)  AU(p,5,T)=q(p)^0{i,T)  +  ^  MklL  (i  -S-p-L.)  q(P> 

n=l  ^  y 


i.,,  (ir  )  r  .  .  2N+3 

whe  re  Rk  <  —N+l -  \  j  Q{\  )  |  \  ;  0  <  p  <  oo  ; 

“  (N+l)  !  *>=0 

and  where  i|/  (§,t  )  is  defined  by 


(C-9)  )  =  \  (T  ')n 


M  2  > 


First;  it  is  noted  that  the  function;  4jq(£;T  ) >  is  the  temperature  response 
for  the  flat  plate  with  a  uniform  input  over  the  entire  front  surface; 
except  for  a  multiplicative  constant: 


(C-10)  T  )  =  T  +  ^2  y 


1  ('- 


2  2  -T ■ 

-  n  TT  T 


COS  nTT  4  -n2  TT  2  T 

- ; -  e 


The  summation; 


;  is  given  by  Jolley.  13  It  should  be  noted 


that  the  integrals  in  Eq.  (C-9)  can  be  expressed  in  terms  of  incomplete 
gamma  functions. 
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Thus  the  heading  term  in  the  small-time  expansion  is  the  product 
of  nhe  local  heat  flux  input  multiplied  by  the  response  for  the  flat  plate 
wi'h  3.  uniform  input  over  the  entire  surface.  This  particular  approxi¬ 
mation  is  val:d  on  the  surface  as  well  as  interior  to  the  slab.  Thus  it 
is  concluded  that  for  time  sufficiently  small  and  the  input  flux  density 
sufficiently  smooth  (as  prescribed  by  conditions  (a)  and  (b)),  the  radial 
diffusion  is  negligible.  A  similar  result  will  be  presented  in  the  next 
section  with  the  smoothness  requirement  weakened.  However,  the  flat 
plate  -  uniform  input  response  approximation  will  hold  only  on  the  input 
surface  and  at  the  center  of  the  front  surface. 


D.  The  Green's  Function 


In  this  section  the  Green:s  function  will  be  examined.  Several 
expansions  will  be  presented,  some  of  which  converge  well  for  t  small 
and  some  cf  which  converge  well  for  r  large. 

Some  additional  general  asymptotic  results  for  the  temperature 
response  as  t  0  are  given.  The  requirements  on  the  input  flux 
density,  q£p),  are  not  as  strong  as  they  were  in  Section  C. 

To  begin,  the  integration  of  Eq.  (B-8)  with  respect  to  \  can  be 
performed.  The  result  is  given  by 


where  I0(x)  the  modified  Bessel  function  of  the  first  kind.  This 
expression  may  now  be  put  in  a  more  compact  and  in  some  respects 
a  more  useful  form  by  introducing  the  integral  form  for  the  Bessel 
func-ion:*- 


(D-2) 


r?  P  P* 

I  •’*’ 
0=  0 


d9 


i  07  .’4 
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Using  this  representation  one  arrives  at  the  following  expression; 


TT  T 


P-3)  G(p 


r. 


0=0  T  ’=0 


0  *  /  —  ,  iir  t  1 


exp 


~  IP  "  P '  2 
4t  1 


where 


fp  -  p  1  2  =  p  2  +  p  12  *  p  1  cos  0 


Making  use  of  Watson*s  suggestion  for  studying  integrals  containing  a 
product  of  two  Bessel  functions  of  the  same  order,  18  as  in  Eq.  (B-8), 
it  is  seen  that  17 


if 

P-4)  JQ(Xp)  J0(\p')  =  ^  f  J0(\|p--p'|) 

0=0 

where  the  quantity,  |p  -  p'|,  is  defined  in  Eq.  (D-3),  above. 

Using  this  representation  in  Eq.  (B-8)  the  expression  for  the 
Green,s  function  becomes 

(D-5)  G(  p,  p«,  £,r)  = 

TT  T  00 

S  $ 

9=0  T  ’=0  X  =  0 

If  one  now  performs  the  X -integration  the  result  is  identical  to  Eq. 
(D-3).  18  The  next  logical  attempt  to  simplify  any  of  the  above  three 
representations  would  be  to  perform  the  -^-integration.  However, 
before  this  a  short  comment  on  the  two  expansions  of  the  theta  function 
(Eq.  (B-6a)  or  (B-6b))  is  necessary. 

If  the  "long-time”  expansion,  Eq.  (B-6a),  is  used  in  either  Eq, 
(D-l)  or  Eq.  (D-3),  or  if  the  short-time  expansion,  Eq.  (B-6b),  is 
used  in  Eq.  (D-  5),  then  the  resulting  integrals  are  of  the  form 


4 


ITT  T 


P  -  P 


-X2  T 
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(D-6) 


It  has  rot  been  possible  tc  express  these  integrals  in  any  simple 
fashion.  The  factor 


has  the  form  of  a  generating  function  for  Bessel  functions,  but  as  yet 
it  has  not  been  possible  to  make  use  of  this  fact. 

Therefore,  the  expansion  of  the  theta  function  which  is  to  be  used 
in  each  of  the  above  equations  is  determined  by  this  restriction. 

If  the  short-time  expansion  of  the  theta  function  is  substituted  in 
Eq.  (D-l)  tho  term-by-term  integration  can  be  performed  only  if  the 
Bessel  function  is  also  expanded  in  its  series.  The  resulting  expression 
is  very  clumsy,  and  because  of  this  complexity  it  has  not  been  possible 
to  determine  the  range  of  convergence  of  the  series.  It  simplifies  for 
special  cases,  but  these  special  cases  can  be  obtained  by  more  direct 
methods , 


The  result  is  given  by  the  expression: 


(D- 7}  G{p.  p',4»t 


1 


t! 


Z  tF~  Ls  (m  !  ] 
m=0 
n= 0 

e  -  +  1 


P  P 


>2  -tp'2  +4  (n+  i  i/Z)z 


Zm 


r  Zm  +  2  , 


t  P2  +P,Z  +4(n+  6  i/Z)z 


4t 


where  F(a 


xfp2  +  P*a  +4(n+-eg/2)i 
00 

,  x)  ~  ^  ta  ^  is  the  incomplete  gamma  function.  19 
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For  the  special  case  p  =  £  =  0,  which  is  of  interest  in  the  case  of 
focused  inputs  as  mentioned  in  the  introduction,  the  above  expression 
simplifies  to  the  following: 


(D-8) 


N 


wl 


n=l 


1 

1 

r(i 

p' 

’  W 

1 

—  r  I 

q  p'2  +  4n2  \ 

^  p2  +  4n2 

1 

v2  4t  / 

+  R3  where 


T r  \|  ITT  _ i  /_ 

R  3  <  — ^ -  e  1/T  , 


N  =  0 


—  _{N+1)2/t  1 


N 


,  N  >  1 


The  result  in  Eq.  (D-8)  also  follows  from  Eq.  (D-3),  by  substituting 
the  theta  function  from  Eq.  (B-6b).  The  gamma  function,  T(j9  x),  is 
related  to  the  incomplete  error  function.  The  relationship  is  given  by 


(D-9) 


where 


erfc  ( 


nTV  erfc  (fj) 

00 

—  f  e't2 

Npr  J 

t  = 


It  will  be  recalled  that  the  incomplete  error  function  enters  into  the 
representation  for  the  small-time  solution  of  the  flat-plate  response 
for  the  case  of  uniform  input  over  the  entire  surface. 

The  flat-plate  response  with  a  uniform  input  over  the  entire 
surface  can  be  deduced  from  the  above  expression  by  letting  a  ~*oo 
while  Q/ita2  is  held  constant. 
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An  asymptotic  expansion  to  one  term  in  terms  of  f T  as  t  -*■  0 
can  be  obtained  from  Eq.  (D-8).  The  result  again  displays  the  local 
nature  of  the  response  for  time  sufficiently  small.  From  Eq.  (D-8) 
one  has 


(D-10) 


A U(0,  0,  r 


+ 


where 


R2  < 


Q 

7k 


The  integral  in  Eq.  (D-10)  is  0(  \pr  )  as  t  -►0;  thus  for  sufficiently 
small  time  it  will  afford  a  good  approximation.  In  general,  however, 
the  integration  in  Eq.  (D-10)  cannot  be  carried  out  in  closed  form. 

It  is  possible,  however,  to  obtain  an  asymptotic  expansion  of 
AU(0,  0,  t  )  as  t  0  to  one  term  under  the  following  assumptions: 

(a)  q(p  )  ~  q(o)  +  0(p  )  as  p  -0 

(b)  q(p  )  <  M,  a  positive  constant,  for  all  p  . 

Under  these  assumptions,  which  are  the  smoothness  requirement 
for  this  case,  one  arrives  at  the  result 


P-11) 


A U(0,  0,  t  ) 


q(o)  +  0(t  ) 


as  t  0 


This  result  states  that  for  r  sufficiently  small  the  response  at  p  =  £  =  0 
is  given  by  the  local  input  multiplied  by  the  small-time  response  of  the 
plate  to  a  uniform  input  over  the  entire  front  surface.  This  is  the 
result  one  would  expect  on  a  physical  basis.  This  result  is  almost  the 
same  as  that  deduced  in  the  previous  section  in  Eq.  (C-5).  And  for 
p  0  and  t  sufficiently  small,  the  one-term  expansion  of  Eq.  (C-8) 
is  the  same  as  (D-ll). 

Next,  attention  is  focused  on  the  integration  of  Eq.  (D-3).  For 
this  expression  the  small-time  expansion  of  the  theta  function  is  again 
used.  The  term-by-term  integration  yields  the  following  expression: 
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+  r3  . 


N  0 

0  <  4  <  1  . 
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This  result  is  quite  similar  to  that  in  Eq.  (D-7),  except  that  the 
summation  on  MmM  in  Eq.  (D-7)  has  been  replaced  by  an  integration 
on  Gin  Eq.  (D-12),  and  it  has  been  possible  to  deduce  a  remainder  for 
this  last  expansion.  This  approach  then  provides  a  method  for  a 
partial  summing  of  Eq.  (D-7).  Furthermore,  using  Eq.  (D-12),  with 
time  sufficiently  small  the  results  in  Eq,  (D-10)  can  be  generalized  to 
yield 


(D~13) 


AU(p  ,  %,  t  ) 


q(p  ')  P’ 


P  ,=o 


f 

0=0 


\l 


p 


V  2  +  e2 


+  R3 


where  the  remainder  is  bounded  for  all  p  by 


R,  < 


Q 


12 


Q 

R,  <  "TT 
3  £k  8 


•1/t 


4=0 


1  -(1-4/2)2/t  -(1  +  £/2)2/t 

—  e  +  e 

3 


The  result  in  Eq.  (D-13)  reduces  to  Eq,  (D-10)  for  p  =  0  since  the 
integrand  becomes  independent  of  Band  the  0-integration  can  be 
performed.  In  this  case  the  dependence  on  £  can  be  retained,  and  in 
this  sense  it  is  more  general  than  Eq.  (D-10).  An  asymptotic 
expansion  can  be  deduced  similar  to  that  given  in  Eq.  (D-ll).  Equation 
(D-13)  is  written  for  p  =  0  as  follows: 


(D-14) 


AU(0,  4,  t 


1 


q(p') 


p'=o 


^  42  +  P 


I  2 


+  R, 


where  R3  is  bounded  as  above  in  Eq.  (D-13). 
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The  asymptotic  development  of  Eq.  (D-14)  yields  (under  the  same 
assumptions  as  those  for  Eq.  (D-ll))  the  result: 


(D-15) 


AU(0,£,r)~q(o)  2 


T 

e-£!/4r.  e  n 

1 - 

1 

V  TT 

2  fr  V2 

'  4  r) 

+  0(t) 


as  r 


At  the  front  surface,  |  =  0,  this  reduces  to  the  result  given  in 
Eq.  (D-ll).  Again  the  local  nature  of  the  response  for  small  time  is 
emphasized.  Furthermore,  it  is  seen  that  this  approximation  is 
independent  of  the  thickness  of  the  slab,  i ,  while  the  remainder  is  not. 


The  next  approach  is  to  perform  the  r  1  integration  in  Eq.  (D-5). 
This  yields  the  following  result: 


t r  00 


(D-16)  G(p  ,  p £,t)  =  -  J 


0=0  X=0 


1  -  e 


-X2  r 


+  2 


n=l 


cos  nTT  £ 


\2  +  n2  TT2 


il  -  e 


-(X2  +  n2  TT  2  )  T 


J0(xfp-p'  I) 


It  has  unfortunately  not  been  possible  to  perform  the  following  integration: 


(D-17) 


R.=  2 


2 


X  =  0  n  =1 


cos  n  ir  £  -( X2  +  n2  it  2 

7^77^  * 


Jq(X  I  P  “  P'  |)X 


However,  an  upper  bound  can  be  obtained  on  this  remainder  and  the  result 
is  an  expansion  which  is  accurate  for  t  sufficiently  large.  The  bound, 
valid  for  all  0  <  £  <  1  and  0  <  p  <  oof  0  <  p'  <  oo,  is  given  by 


( D-18) 


R4  < 


1_ 

6 


T 
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When  the  exponential  factor  is  not  present  it  is  possible  to  perform 
the  \ -integration20  in  the  second  term  of  Eq.  (D-16), 


This  yields  the  following  expansion: 


P-19) 


g(p,p\  e,r)  = 


0=0 


0=0  ^  =  0 


00 


n=l 


cos  mr  £  K0  (nir 


JG(^  Ip  -  p'|)^ 


P'|)  +  R4  . 


It  is  also  possible  to  rewrite  this  last  sum,  making  use  of  the  addition 
theorem,  21  as 

t r 

(D-20)  i  C  K0(mr  |  p  -p'|)  = 

9-0  i0(niTP)  K0Np')  p  <  p' 

Io(mrp')  K0(mrp  )  p  >  p  ' 


It  appears  that  the  only  way  to  handle  the  first  integral  is  to  make 
use  of  the  equality 


T 


By  making  use  of  this,  both  the  X -integration  and  the  r  '-integration  can 
be  carried  out  in  terms  of  the  incomplete  gamma  function,  r(0,x),  or 
the  exponential  integral. 

Incorporating  these  suggestions  and  using  Eq.  (D-20)  in  Eq.  (D-19), 
one  arrives  at  the  following  expression  for  the  Green's  function,  which  is 
useful  for  t  large: 
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(D-22) 


G(p».p  1 , 5,  t) 


2tt 


r 


0=0 


+ 


z 


n=l 


I 

n-1 


cos  mr  £l0(mr  p )  K0(n7rp'  ),  p  <  p1 


cos  n.TT  ^I0(nTT p1  )KQ(n7Tp  ),  p  >  p’ 


1  e 


-TT  2T 


+  R4  ,  where  R4  <  — 


0  <  p  <  oo 
0  <  p'<  oo 

0  <L£  <  1 


It  has  not  been  possible  to  deduce  any  general  asymptotic  results 
for  r  large  similar  to  those  in  Eqs  .  ( D - 11 )  and  D-15)  for  t  small. 

In  this  section  several  expansions  have  been  presented  for  the 
Green's  function.  Some  of  these  converge  quickly  for  t  small,  and  some 
of  them  for  r  large.  In  addition,  some  general  asymptotic  results  for 
the  temperature  response  as  t  -►  0  have  been  given.  The  requirements 
on  q(p)  needed  to  obtain  these  results  are  not  as  strong  as  in  the  case  of 
the  development  in  Section  C. 


E *  Gaussian  Input 

The  thermal  response  to  the  Gaussian  input  is  discussed  in  this 
section.  Two  small-time  expansions  for  the  temperature  at  the  center 
of  the  spot,  p  =  £  =  0,  are  given  and  their  relative  accuracy  is  compared. 

In  addition  the  criteria  necessary  to  use  the  asymptotic  expansion  pre¬ 
sented  in  Eq.  (C-8)  are  satisfied  for  this  input  flux  density.  For  the 
parameter  4r/r\^  <  1,  this  yields  a  convergent  series  which  can  be 
summed  to  give  a  form  convenient  for  numerical  computation.  This 
result  can  also  be  obtained  by  other  techniques.  The  asymptotic  expansion 
deduced  from  Eq.  (C-8)  is  of  course  not  restricted  to  the  spot  p  =  £  =  0. 
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The  Gaussian  input  is  defined  as  follows: 


(E- 1)  q(p )  -  £fk  ‘--z  ■  exp(-p2  /no2  )  ,  p>  o 

1  O 


where  Q  is  the  total  heat  flux  input  and  r\  =  &/£  where  naM  is  the  ’’radius" 
of  the  spot.  The  transformed  input,  Q(^;  is  given  by  a  function  of  the 
same  form ,  2  2 


(E-2)  Q(M~  Qj-  ■  exp>(-\2  rf  /4). 

2  7T  £  K  ° 


The  fact  that  the  parameter  X  enters  the  expression  for  the  temperature 
field  in  the  same  manner  as  the  factor  exp(-X2  t1)  in  Eq.  (  B-5)  makes  this 
input  particularly  amenable  to  analysis.  Here  is  the  first  divergence  from 
the  rather  general  results  obtained  in  the  last,  two  sections.  It  is  a  fruitful 
divergence  because  of  the  special  form  of  the  heat  flux  input.  The  short- 
time  expansion  of  the  theta  function  given  in  Eq.  (B-6b)  is  going  to  be  used 
with  the  F{p ,  t;  X)  form.  This  will,  of  course,  limit  the  results  to  the 
case  £  ~  0  and  a  one-term  expansion.  This  is  because  of  the  difficulties 
presented  by  integrals  of  the  form 


[aT^b  ) 

g(r')  e  T  : 


0 


as  was  discussed  above  in  Section  D.  An  upper  bound  on  the  remauuig 
terms  in  determined.  The  result  is  given  by  the  expression 


(E-3) 


ALT(  0 , 0 ,  t  ) r 


Q_ 

ak 


Rt 


whe  re 
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_Q_  f  4kt/az  A  -1/t 

Rs  <  fk  12  V 1  +  4kt^a  y  6 


There  are  several  things  to  note  about  this  expression.  First, 
for  r  small,  this  approximation  is  best,  that  is  the  error  bound  is 
controlled  by  the  term  e"^/T.  Second,  the  approximation  itself  does 
not  depend  on  the  thickness  of  the  slab,  but  error,  R5  ,  does.  Third, 
for  t  fixed  the  approximation  becomes  better  as  the  thickness  of  the 
slab  increases.  This  statement  is  somewhat  redundant  since  this  con¬ 
clusion  follows  from  the  first.  However,  the  conclusion  in  this  form 
emphasizes  the  nsmall-timen  aspect  of  the  approximation.  Furthermore, 
as  l  *-►  oo  this  response  approaches  that  of  the  semi-infinite  body  at 
p  =  £  =  C  since  the  remainder  R5  0  as  I  oo  .  A  plot  of  this  approxi¬ 
mation  is  shown  in  Fig.  2;  while  the  plot  appears  to  say  that  the  tempera¬ 
ture  is  bounded  as  t  -*•  oo  ,  this  is  not  the  case  since,  as  the  remainder 
indicates  the  approximation  eventually  becomes  very  poor  as  t  increases 
with  the  other  parameters  fixed. 


Fig.  2.  The  response  to  the  Gaussian  input  as 
given  by  the  arc-tan  representation. 
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The  second  expansion,  deduced  from  the  large-time  theta  function 
expansion  given  in  Eq,  (B-6a)  and  the  F(p,  £,  tt  ;X  )  form,  is  convergent  for 
all  t  .  In  contrast  to  the  arc-tan  representation,  the  first  term  of  this 
expansion  does  depend  on  the  thickness  of  the  slab,  ft  .  The  expression 
for  the  remainder  shows  that  the  remainder  is  again  small  for  4kt/a2 
small,  but  by  contrast  to  the  arc -tan  representation,  the  representation 
is  best  when  the  parameter  r)0  is  large.  The  representation  is  given  by 


(E-4) 


*J(0,  £,t)  = 


Q 


Ztt  l k 


In 


J1 


+  4Kt 


+ 


N 


Z 

n=l 


cos  nTT  £ 


„  2  _  2  „2 
n  u 


)  -r  (°- 


_  2  2  2 
n  TT  p 


1  + 


•  exp 


2  2  2 
n  tt  6  t|q 


4 


+  ^6 


where 


From  the  above  expression  for  the  remainder  it  can  be  shown  that 
for  all  t  the  remainder  can  be  made  arbitrarily  small  by  increasing  the 
number  of  terms,  N,  in  the  expansion.  However,  as  N  becomes  large  the 
remainder  vanishes  only  as  1/N  and  the  convergence  is  very  slow.  However, 
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in  contrast  to  the  arc-tan  representation,  arbitrairly  good  accuracy  can 
in  principle  be  obtained  from  this  expansion  for  all  t  .  Finally,  the  accuracy 
of  this  representation  when  only  one  term  is  retained  is  not  as  good  as 
that  obtainable  from  the  arc-tan  representation,  as  an  examination  of  the 
remainders  show.  A  graph  of  the  logrithm  representation  is  given  in 
Fig.  for  two  values  of  rj  ,  using  the  log  term  plus  the  first  term  in  the 
s  e  ries . 


0.002  0.004  0.006  0.008  0.010  0.012  0.014 


4  k  t 
02 


Fig.  3.  This  is  the  response  to  the  Gaussian  input  as  given 
by  the  log  representation  using  the  log  term  plus  the 
first  term  in  the  series.  The  response  is  plotted  for 
two  values  of  t| 

This  graph  shows  that  for  t  and  the  total  heat  input  fixed,  the 
temperature  decreases  as  p  increases.  This  is  the  expected  result. 

The  requirements  necessary  to  use  the  asymptotic  expansion  in 
Eq.  (C-8)  are  fulfilled  by  the  Gaussian  input.  The  first  three  terms  in 
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the  expansion  are  given  by 


AU  (P;  T  )  = 


q(P) 


+0(^i  T  )  + 


where  the  remainder  is  bounded  by 


t  )> 


for  allp.  This  may  also  be  written  for  the  case  of  the  Gaussian  input  as 
follows . 1 2 


(E-5)  AU(p,e;  t) 


whe  re 


(£  A  T  ) 
n 


J  (t  *  )  63  (  ji-  ,  i  IT  T  '  )  , 

t'  =  0 


Ln  (x) 


Laguerre  polynomials! 


2  3 


and 
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Rb  ^  q(p ) 


N+l 


JN+1 


+N+1  (£»  T)>  N  >  0. 


The  expression  for  the  upper  bound  can  be  simplified  by  use  of  the  fact 


that 

(E-6) 

+N+1  (£»  t)  <  tN+1 

and 

Ln/£^<  ep2/< 

where  (  ^ ,  t )  is  given  by  Eq.  (C-10). 

One  last  representation  for  the  temperature  fields  will  be  given. 

It  can  be  obtained  by  a  variety  of  techniques.  One  technique  is  by  in¬ 
troducing  the  integral  definition  of  ^n(£>T)  into  Eq.  (E-5)  and  then  noting 
that24 


00 

(E-7)  y  znL  (x)=-i —  exp  , 

Li  1  -z  z  -1 


Z  <  1  . 


n  =  0 


Another  method  for  arriving  at  the  same  result  is  by  direct  integration 
of  Eq.  (B  -5) : 


(E-8) 


aj(p.e.T) 


\2 

J0(Xp)K  . 


The  integration  with  respect  to  X  can  be  performed17  to  yield  the  following 
result : 

Q  fT 

(E-9)  AU(p,e.T)  =  77^  J  03 

r'  =  0 


—  ,  i  TT  T 
2 


4t  1  + 


exp 


4t  1  + 


2 

^o 
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This  form  of  the  response  seems  most  useful  if  a  numerical  computation 
is  required.  The  term  q(p)  can  be  factored  from  this  expression  to 
yield  the  result  indicated  by  the  sum  in  Eq.  (E-7). 


F .  The  Uniform  Input 


In  this  section  the  temperature  fields  due  to  a  uniform  input  on 
a  circular  spot  of  finite  radius  are  presented.  The  input  is  defined  in 
Eq.  (F-*l)  below.  Because  of  the  discontinuity  in  q(p).  the  results  for 
this  case  are  not  as  abundant  as  they  were  for  the  case  of  the  Gaussian 
input* 


One  of  the  results  is  an  extension  of  Oosterkamp's  work.3  It  is 
extended  to  include  results  for  0  <  e<  l  and  p  =  0,  while  his  results 
held  only  for  p  =  £  =  0.  Included  is  an  upper  bound  on  the  remainder 
which  serves  to  give  a  precise  meaning  to  the  approximation.  Graphical 
information  is  also  included. 

The  long-time  Green's  function,  defined  in  Eq.  (D-22),  is  used  to 
obtain  the  response  for  the  uniform  input.  A  simple  explanation  in  terms 
of  the  effect  of  radial  diffusion  is  given.  Graphical  information  is  also 
included. 

The  results  of  this  section  also  show  that  when  the  spot  size  in¬ 
creases  to  infinity  and  the  input  flux  density  is  held  constant,  the  tem¬ 
perature  response  approaches  that  of  the  flat  plate  with  a  uniform  input 
over  the  entire  front  surface.  This  is  of  course  to  be  expected. 

The  uniform  input  is  defined  by  the  expression 


(F-l) 


q(p)  = 


i 

k  ’ 


p  <  % 


=  o 


p  > 
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where  a  is  the  radius  of  the  spot  and 


The  transformed  input  is  given  by  the  expression25 


(F-2) 


Q(x  ) 


l  Ji(M0) 

k  '  T1°  X 


The  Green‘s  function,  G(p ,  p* ,  t)  seems  to  be  the  most  useful  for  this 
input  because  the  integrand  of  the  expression, 


(F-3) 


4U(p.  it  t) 


J  Q(X)  F(p,  £9t;\)\  , 
\  =0 


contains  the  product  of  two  Bessel  functions  of  different  order  and  argument, 
and  these,  as  mentioned  before,  are  the  most  difficult  to  handle.  The 
only  possible  approach  using  the  representation  in  Eq.  (F-3)  is  to  make 
use  of  the  expansion2-6 


(F-4) 


Ji  tiQ)J0(Xp)  = 


00 

Mo  y  (-Dm 

2  Z  (m!)2 

m-0 


P  >  po 


Xtio  y  <-i)m 
2  L.  m !  (m+1) ! 
m  =  0 

P  <  V 


Mo 


cm 


-rn,  -l-m;l;p2  / r| 


2^ j  (  -m,  -m;  2;  r|0/p2  ) 
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The  functions  2  F2  (-m,  a;b;x2  )  are  hypergesnetric  functions,  27  ’  and  be¬ 
cause  of  the  factor  (-m)  are  polynomials  of  degree  in  in  x2  .  This  ex¬ 
pansion  allows  one  to  perform  the  \  integration,  but  the  integration  on 
t*  cannot  be  carried  out  for  £  =  0.  The  special  case  p  =  0  can  be 
handled  with  the  form  given  in  Eq.  (F -  3)  but  will  be  obtained  later  with 

This  result  is  given  by  the  expression 

L-f  <i-.-,'o/4T')e, 

T  *=  0 

Before  passing  to  the  Green* s  function  approach,  several  comments  about 
this  result:are  in  order.  First,  it  is  recalled  that  the  solution  to  the  flat- 
plate  problem  with  a  uniform  input  over  the  entire  surface  is  given  by 

(F-6)  AU  =  — — 

ir  az  k 


T  *  =0 


the  use  of  the  Green's  function. 


(F-5)  A U(0,  £,t)  = 


Q 


t t  a2 


This  is  precisely  the  result  that  one  would  obtain  from  Eq.  (F-5) 
by  holding  Q/tt  a2  constant  and  letting  a  -*oo.  Also  for  t  small  the 
effect  of  the  exponential  term  in  the  integra]  of  Eq.  (F-5)  is  small,  and 
thus  it  is  seen  that  when  t  is  small,  the  response  of  the  temperature 
at  p  e  0  is  that  of  the  flat  plate  with  a  uniform  input  over  the  entire 
surface.  This  again  indicates  the  local  nature  of  the  response  for  t 
small. 


If  the  first  term  in  the  small-time  expansion  of  the  theta  function 
for  £  =  0  is  retained  and  the  integration  performed,  the  result  is  the  same 
as  Oosterkamp* s.  It  is  not  necessary,  however,  to  require  £=  0.  The 
integration  can  be  performed  even  for  £  £  0.  Also  an  upper  bound  on  the 
remainder  can  be  computed.  For  £  =  0  the  result  is  given  by 


(F-7) 


AU(0,  0,  t)  = 


_Q _ 

tt  ak 


1  4T 


2 


r 


+  r2 
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where 


R,<-“  ,E  e-l/T 
2  £k  12 


Attention  is  called  to  the  similarity  of  this  result  and  the  small-time 
approximation  given  in  Eq.  (D-15). 

The  result  in  Eq.  (F-17)  may  also  be  rewritten  as  follows: 


(F-8) 


AU(0,  0,  t)  =  -9_ 

it  ak 


4Kt 


,  a2  /  4Kt 
1  -  e _ 


+ 


a 

2  Kt 


+  R2  » 


which  is  closer  to  Oosterkamp's  notation.  A  plot  of  this  result  is  shown 
in  Fig.  4.  These  formulas  emphasize  the  "small-time"  aspect  of  the 
approximation.  Note  that  the  approximation  does  not  depend  on  the  thick¬ 
ness  of  the  slab  (it  is  therefore  the  response  of  the  s emi -infinite  slab 
since  R2  -*■  0  as  i  -*oo),  but  the  bound  on  the  remainder  does  depend  on 
the  slab  thickness.  The  requirement  that  Eq.  (F-7)  or  (F-8)  be  a  good 
approximation  is  that 


T 


Kt 


be  small. 


Thus  for  any  fixed  t  (in  seconds),  the  thickness  of  the  slab  can  be  in¬ 
creased  so  that  the  approximation  in  Eq.  (F-7)  or  (F-8)  offers  an  approxi¬ 
mation  as  good  as  needed.  Some  typical  values  of  the  normalized  remainder 
R/[Q/fk],  are  given  in  Table  I  below. 

TABLE  I 

Some  typical  values  of  the  remainder 
associated  with  Eqs ,  (F-7)  and  (F-8). 


T 

R/l 

[Q/fk] 

.  01 

5  . 

10'46 

0.1 

2  • 

O 

1 

1.0 

5  • 

1CT2 
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4  kt 
a2 


Fig.  4.  This  is  the  response  to  the  uniform 
stop  as  obtained  by  Oosterkamp. 

For  0  the  response  is  given  by  the  following,  which  is  quite  similar 
to  the  earlier  result  in  Eq.  (F-8)  : 


whe  re 
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£  <  Q_  Jttt 

7  tk  ’  4 


o  <  £  <  l. 

Note  again  that  there  is  no  dependence  on  the  thickness  except  in  the  re¬ 
mainder,  and  that  as  I  -^oo,  R?  0 .  For  £  =  0  the  above  approximation 
reduces  to  the  former  case. 

Now  the  long-time  Green’s  function  is  used  to  obtain  the  temperature 
response  for  t  large.  This  Green's  function  is  given  in  Eq.  (D-22).  In 
order  to  obtain  this  result,  use  is  made  of  the  two  integrals  :28>  29 


-O-S/2)2 


+  i 


-(i+e/2)2 


rp 

\  P  'I0(nirp  ')  =  -£-  I  (mrp), 
J  nrr 


p'  =  0 


\  p  'K0(nrrp  ')  =  _L  [pK,  (mrp)  -  T|  Kt  (nirr)  )  ] 
ntr 


P  ,=P 


The  temperature  fields  for  p  <  t|  are  given  by 


(F  -10) 


4U(p,£,t)  =  L  J_  f f  P'  r( o, 

TTa  k  Ztt  J  J  y 


P  P  '  |  2 

~4t 


p’  =  0  0  =  0 


,  Q  f  V7  cos  niT  £  r 

+  2  - 2  -  >  -  [K0(mrp)pl  (mTp( 

na  k  ^  nir  1 

n=l 


+  I0(nirp)  {pK,  (ntrp)-r|0K1  (nTrr|0)}  ]+Ri  . 
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The  temperature  fields  for  p  >  rj  are  given  by 


(F-ll)  OJ(P,e,T)=  -9-y  i  J- 

ira  k  2 


•  r.f 

p 1 =o  0=0 


r  o. 


p-p 


'  2 


4t 


where 


+  2 


oo 

Q  i  cosrnr  g 


it  a  k 


Zcosmr  5 

- 'K0(nnp-dr)oIi  (nirr|o)  +  Rj 


n=l 


and 


p-p1  2  =  p2  +  pi2 


2p  p  1  cos  0 


R  <  A  ® 


i  fk  lZn- 


Some  typical  values  of  the  normalized  remainder  are  given  in  Table  II. 


TABLE  II 


Some  typical  values 

of  the  normalized 

remainder  for  Eqs . 

(F-10)  and  (F-ll). 

T 

R  /[Q/lk] 

0.1 

.*10 

1.0 

.4  •  10"3 

10. 

.  15  •  10-5 

Again  attention  is  called  to  the  fact  that  as  a  -►  oo  but  with  Q/ira2 
held  constant,  the  above  result  (Eq.  (F-10)  must  be  used  because  of  the 
restriction  p  <  tj  )  reduces  to  the  response  for  the  flat  plate  with  a  uniform 
input  over  the  entire  front  surface,  except  for  the  term  which  contains 
both  a  £  dependence  and  the  time -dependence .  This  term  has  been 
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absorbed  by  the  remainder.  The  last  summation  becomes  independent  of 
p  since  it  is  p  multiplied  by  the  Wronskian  of  IQ(x),  KQ(x).  The  initial 
integration  can  be  performed  and  is  proportional  to  t  but  is  also  in¬ 
dependent  of  p  . 

Again  the  special  case  p  =  0  can  be  given  more  attention.  Setting 
p  =  0  in  Eq.  (F-10)  one  arrives  at  the  following  result: 


(F-1Z) 


AU( 0,£,  r ) 


+ 


£  + 


3 


;  l 

n=l 


COS  HTT  £ 


K,  (nim0) 


+  R, 


where 


Again  as  ru  -*oo  but  - 2  remains  fixed,  this  result  reduces  to  the  flat- 

&  'o  it  a 

plate  response  with  a  uniform  heat  flux  density  input  over  the  entire  surface. 


Note  that  the  first  term  in  this  response  is  mtm  multiplied  by  a 
factor  in  brackets.  This  factor 


§ 


+ 


i  rfo, 

4t  ^  4t  J  J  , 


will  be  called  the  diffusion  factor,  and  it  carries  the  effect  of  not  having 
a  uniform  input  over  the  entire  surface.  Note  that  it  is  independent  of 
the  slab  thickness.  This  factor  is  shown  graphically  in  Fig.  5.  This 


1107-14 


35 


Fig.  5.  This  is  a  plot  of  the  diffusion  factor  defined 
in  Eq.  (F-12)  assoicated  with  the  long-time 
response  for  the  uniform  spot. 


figures  shows  that  as  naM  increases^  for  t  fixed^  the  effect  of  radial 
diffusion  at  p  *  0  decreases*  Also  it  shows  that  as  t  or  K  increases 
for  MaM  fixedj  the  effect  of  radial  diffusion  at  p  *  0  also  increases* 
The  last  term  in  the  response  of  Eq.  (F-12)> 


n  *  1 


9 


is  independent  of  time^  and  as  t]  ***  00  the  last  sum  vanishes  and  also  the 
response.  Thus  the  larger  rj0^  the  less  the  effect  of  radial  diffusion  at 
p  =  0.  Fig.  6  displays  this  last  sum  for  a  range  of  r|Q. 
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Fig.  6.  The  last  summation  in  Eq.  (F-12)  is 
shown  for  £  *  0  and  0.  1  <  tt  <60. 


G.  Conclusions 


Transformation  techniques  are  used  to  obtain  several  alternative 
representations  of  the  solution  in  terms  of  the  Green*s  function  and  its 
transform.  Some  representations  converge  well  for  time  small,  and 
some  for  time  large. 

Several  asymptotic  expansions  are  deduced  for  the  small-time  response 
of  the  temperature  fields  under  various  assumptions  on  the  heat  flux  input 
density. 


An  extension  of  the  results  of  Oosterkamp  has  been  derived  as 
well  as  a  verification  of  the  validity  of  the  approximations  through  the 
derivation  of  an  upper  bound  on  the  error. 

Two  specific  inputs  have  been  considered  in  detail:  the  Gaussian^ 
and  the  uniform  spot.  These  two  inputs  can  serve  as  models  of  the  heat 
flux  input  caused  by  such  sources  as  the  laser;  the  study  of  the  thermal 
response  is  the  initial  step  in  a  study  of  the  thermal-elastic  response 
of  a  slab  to  an  input  such  as  the  laser. 
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